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ABSTRACT 

A relatively  simple  double  pass  method  is  used  to  (a)  reduce 
longitudinal  spread  of  satellite  navigation  fixes  and  (b)  obtain  antenna 
height  or  sea  level  height  on  board  a ship  in  a fixed  harbor  or  fixed  open- 
sea  position  (over  an  acoustic  beacon).  Required  are  the  satellite 
navigator  output  of  latitude,  longitude,  elevation  angle,  and  path  geometry 
for  two  consecutive  passes  (one  E and  one  W of  the  site)  of  the  same 
satellite . 

The  standard  deviation  of  one  such  double  pass  observation  in  a larger 
data  set  is  usually  below  20  m in  both  longitude  and  height.  Improvements 
in  the  position  accuracy  of  a site  (as  measured  by  a , the  standard 
deviation  of  the  mean  or  the  957,  confidence  interval)  are  obtained  by 
restricting  observations  to  passes:  with  a high  number  of  symmetric  Doppler 

counts,  between  given  elevation  angles,  with  a small  number  of  iterations  of 
the  satellite  navigator,  and  with  eliminations  of  double  pass  results  that 
are  outside  specified  limits  (usually  2ci)  in  height  and/or  longitude.  An 
accuracy  of  10  m or  less  is  usually  achieved  within  four  to  five  days  of 
recordings . 

Sea  level  height  thus  determined  on  Board  the  R.  V.  Kana  Keoki  in  harbors 
around  the  Pacific  and  on  board  the  Glomar  Challenger  in  position  over  Deep 
Sea  Drilling  Project  sites  in  the  Atlantic  show  good  agreement  with  geoidal 
height  models. 

Determining  sea  level  (with  respect  to  the  Navigation  Satellite  Reference) 
using  the  double  pass  method  thus  provides  an  important  independent  alternative 
to  Sky  Lab  Radar  Altimetry  for  the  measurement  of  geoidal  height  over  the  open 
oceans. 
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The  Navy  Navigation  Satellite  System  became  operational  in  1964. 

Since  then,  commercial  equipment  has  become  available  with  which  position 
fixes  in  longitude  and  latitude  can  be  obtained.  Descriptions  of  the 
satellite  system  have  been  given  by  Stansell  (1970,  1976).  In  the  1976 
description  an  error  budget  is  presented  (on  p.  481)  for  different  sources 
of  uncertainties  at  a fixed  location,  and  the  "root  sum  square"  of  the 
errors  is  given  to  be  in  the  range  of  18  to  35  meters  for  "laboratory  type" 
measurements,  whereas  field  results  usually  lie  in  the  range  of  27  to  37  m. 
These  values  apply,  according  to  Stansell,  after  the  time  when  the  polar 
motion  error  was  modeled  and  included  as  an  adjustment  to  the  transmitted 
orbit  parameters.  The  errors  found  during  the  present  study  (standard 
deviation,  a ) are  about  20  m or  less  in  both  height  and  longitude 
when  results  of  two  consecutive  passes  of  the  same  satellite  in  the  same 
direction  (N-  or  S-bound)  are  combined  into  a simple  "double  pass"  solution 
obtained  from  a fixed  position  on  board  the  R.  V.  Kana  Keoki  when  the  vessel 
was  tied  to  a pier  in  various  Pacific  harbors.  Berg  (1975)  has  shown  in  a 
detailed  example  that  the  longitude  positions  obtained  by  the  double  pass 
method  are  more  accurate  than  those  obtained  by  simple  arithmetic  averaging, 
as  was  done  by  Daugherty  (1974) . The  simple  arithmetic  means  show  increased 
standard  deviation  with  increased  error  in  initialized  antenna  height,  an 
effect  that  is  eliminated  in  the  longitude  determination  by  the  double  pass 
method.  Some  very  limited  data  for  fixed  land  stations  in  Colombia  during 
spring  1973,  also  analyzed  by  the  double  pass  method  (Berg,  1975,  and  in 
press),  seem  to  indicate  that  the  standard  error  is  of  the  order  of  9 m, 
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wlthout  application  of  corrections  for  tropospheric  refraction.  These  data 
also  were  obtained  with  standard  commercial  equipment  (a  Magnavox  706). 

The  writer's  interest  in  the  problem  of  SAT  NAV  fixes,  and  especially 
the  antenna  height  correction  problem,  was  generated  through  his  participa- 
tion in  Project  Narino,  a seismic  crustal  study  in  Colombia  that  was  depend- 
ing upon  satellite-defined  positions  reported  by  Woollard  and  Thompson 
(1974)  for  the  location  of  land  shot  points  and  recording  sites  that  would 
be  compatible  with  offshore  shot  points  obtained  by  SAT  NAV  on  board  the 
R.  V.  Kana  Keoki  (Berg,  1975,  and  in  press). 

For  the  present  study,  computer  programs  have  been  developed  with  the 
aid  of  Duncan  M.  Chesley  to  replace  the  graphical  solutions  and  hand 
calculations  that  were  used  to  prepare  the  report  of  Berg  (1975).  These 
programs  allow  various  selection  criteria  to  be  applied  to  the  raw  data  as 
obtained  directly  from  the  satellite  navigator  output  at  a fixed  location. 

Among  these  selection  criteria  are  the  elevation  angle  of  the  satellite 
at  closest  approach,  the  number  of  iterations  of  the  internal  SAT  NAV  computer 
to  obtain  the  latitude  and  longitude  fix,  and  the  number  of  symmetric  Doppler 
count  intervals.  In  addition  to  these  direct  criteria,  selection  from  a 
given  data  set  can  be  made  by  day  or  satellite  number  and  "erratic"  values 
rejected  from  the  double  pass  results.  The  latter  is  achieved  by  calculating 
the  mean  values  of  height  and  longitude  and  their  standard  deviations, 
eliminating  all  double  passes  outside  a specified  height  limit  (usually  2a 
or  3a),  then  eliminating,  in  addition,  double  passes  for  which 
the  longitude  is  outside  the  2a  or  3a  limit  in  longitude.  For  the 
remaining  passes  the  same  procedure  is  repeated  as  many  times  as  necessary 
until  no  values  are  outside  the  specified  limits.  This  last  procedure  works 
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very  well  for  larger  data  sets  but  sometimes  misses  obviously  "bad"  values 
for  small  sets  (five  to  ten  double  pass  determinations),  even  with  a 2 a 
limit.  Results  of  each  such  interaction  step  are  printed  out  by  the 
computer.  Two  output  plots  are  also  obtained  for  each  site:  (a)  all  double 

pass  solutions  and  (b)  the  retained  solutions  (if  different  from  the  set 
under  (a)). 

The  accuracy  of  the  position  data  in  height  (and  longitude)  will  be 
discussed  on  the  basis  of  the  larger  data  sets.  The  height  determinations 
are  particularly  useful,  since  no  geoidal  height  changes  (of  the  order  of 
the  accuracy  of  the  determination)  are  expected  between  visits  of  the 
Kana  Keoki  to  the  same  Pacific  ports  over  the  years  (especially  Honolulu, 
abbreviated  to  HNL  or  HON,  and  Suva  Fiji).  No  comparison  will  be  made 
between  the  determined  latitudes  and  longitudes  and  the  position  of  the 
ship  (antenna)  on  the  maps  of  the  different  locations,  as  had  been  done  for 
Colombia  (Berg,  1975,  and  in  press),  since  this  is  the  specific  subject  of 

a thesis  by  Valerie  Hanna  of  Hawaii  Institute  of  Geophysics,  who  graciously 

. 

provided  the  IBM  data  cards  with  the  raw  data  (the  individual  SAT  NAV  fixes) 
for  all  Kana  Keoki  data  prepared  for  her  thesis.  She  also  prepared  the  data 
cards  for  SAT  NAV  fixes  of  the  Deep  Sea  Drilling  Project  (DSDP)  sites  for 
use  in  the  present  study. 

The  other  result  of  this  investigation  is  the  determination  of  sea 
level  height  or  geoidal  height  at  the  circumpacif ic  ports  visited  by  the 
Kana  Keoki  and  in  the  open  sea  at  the  drill  sites  of  the  Glomar  Challenger. 

On  the  earlier  DSDP  sites  only  a few  fixes  had  been  taken  and  geoidal 
heights  are  known  there  with  only  limited  accuracy.  During  the  later  part 
of  1974,  however,  a larger  number  of  passes  were  recorded  at  each  site,  and 
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geoidal  heights  are  better  determined.  The  accuracy  varies  from  site  to 
site  (as  measured  by  the  standard  deviation  (a)  of  a single  double  pass 
observation).  This  variation  from  site  to  site  is  due  to  residual  motion  of 
the  ship  over  the  drilling  site  between  satellite  passes.  In  general, 
however,  the  cr's  obtained  are  similar  to  those  found  on  board  Kana  Keoki 
in  harbors,  and  the  standard  deviations  of  the  mean  (o^)  of  the  height 
determination  at  the  open-ocean  sites  is  of  the  order  of  12  m or  less,  with 
some  near  2 and  5 meters . 

This  height  determination  at  open-sea  sites  gives  an  important  "ground 
truth"  measurement  with  accuracy  comparable  to  or  better  than  that  of  the 
(limited)  determinations  obtained  by  the  Sky  Lab  Radar  Altimetry  during  a 
single  orbit  (Leitao  and  McGoogan,  1974;  King-Hele,  1976)  and  seems  to  permit 
a direct  comparison.  Since  many  more  passes  and  different  satellites  are 
involved  in  the  open-sea  determinations  presented  here  as  compared  to  the 
observations  of  Sky  Lab,  errors  of  orbital  parameters  are  considerably 
reduced.  King-Hele' s Figure  9 (1976)  shows  differences  between  the  sea 
level  height  determined  from  the  Sky  Lab  altimeter  and  the  Goddard  Space 
Flight  Center  Earth  Model  8 (GEM  8)  of  the  order  of  20  m,  west  of  Australia 
and  also  southeast  of  the  Izu  trench  (SSE  of  Honshu).  It  appears,  therefore, 
that  the  method  used  here  provides  geoidal  heights  from  existing  data  that 
had  not  been  previously  utilized.  The  method  can  also  be  applied  to  future 
data  sets.  The  data  presented  here  suggest  that  a two -day  position 
recording  from  all  operational  Navy  navigation  satellites  can  yield  a 
height  determination  with  om  of  5 m or  less  on  board  the  Glomar  Challenger 
when  located  in  a fixed  position  at  sea. 


No  height  corrections  for  earth  or  ocean  tides  have  been  included  in 
the  calculations.  Their  inclusion  could  further  reduce  the  a in  harbors 
where  the  ocean  tides  are  relatively  large  (such  as  3 to  5 m along  portions 
of  the  coast  of  South  America) , and  where  a seems  to  be  larger  than  in  HNL 
where  tides  are  usually  only  of  the  order  of  .6  or  .5  meter. 


DOUBLE  PASS  METHOD 

Application  to  longitude  and  height  determination 

In  this  section  the  method  will  be  reviewed.  During  a specific  pass 
of  a navigation  satellite,  the  Doppler  counts  obtained  are  predominantly  a 
function  of  the  position  of  the  receiving  antenna  along  the  line  of 
satellite  motion  and  its  distance  from  the  orbital  plane.  Since  the 
navigation  satellites  are  in  polar  orbits,  the  along-track  position  relates 
to  latitude  and  the  cross-track  distance  relates  to  longitude  and  the 
height  of  the  antenna.  Since  distance  is  one  of  the  factors  determined 
from  the  Doppler  counts,  a priori  knowledge  of  the  antenna  height  is  necessary 
to  determine  the  longitude  fix.  Both  "antenna  height"  and  "satellite  height" 
are  taken  as  the  height  above  or  below  the  reference  ellipsoid  or  satellite 
datura.  Figure  1 illustrates  the  geometry  of  two  satellite  passes  in  the 
vertical,  E-W  oriented  (cross-track)  plane  at  the  closest  approach  of  the 
satellites.  The  figure  indicates  that  erroneous  fixes  and  L2  are 
obtained  if  the  initialized  antenna  height  is  different  from  the  true 
antenna  height.  The  amount  of  the  fix  error  L^-L  or  L2-L  depends  on  the 
initialization  error  dH  and  is  a function  of  the  elevation  angle  k of  the 
satellite  at  closest  approach.  For  the  calculations  in  this  report,  this 
function  of  elevation  angle  k:  f(k)  was  taken  from  the  graph  in  Stansell 
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(1976,  p.  486).  It  was  carefully  measured  a-t  the  10°-intervals  marked  in 
the  graph  and  at  65°  and  75°.  The  ratio  of  fix  error  to  height  error  (in 
Stansell,  1976)  was  converted  to  a fix  error  (in  nautical  miles)  per  m 
antenna  height  error,  which  is  the  same  as  the  fix  error  in  minutes  of  arc 
per  ra  antenna  height  error  at  the  equator  and  is  then  f(k).  A smooth 
curve  through  the  data  points  thus  obtained  yielded  the  other  values  at  1° 
elevation  angle  intervals.  The  nautical  mile  was  taken  to  be  1853.2  m. 

Since  fixes  of  the  satellite  navigator  are  given  in  degrees  and  minutes, 
conversion  of  the  fix  error  (to  ra)  can  be  obtained  by  multiplying  the 
dL£(=L^-L)  by  1853.2  • cos  (LAT) , where  the  units  of  dL  are  taken  in  minutes 
of  arc.  Returning  to  Figure  1,  it  is  seen  that  if  f(k)  is  the  fix  error 
(in  nautical  miles)  resulting  from  a 1-raeter  height  error,  then  a dH-meter 
height  error  will  result  in  dH  • f(k)  nautical  mile  error  in  longitude  or  an 
error  of  dH  • f (k)/cos(IAT)  in  units  of  minutes  of  arc. 

If  a number  of  satellite  passes  are  recorded  at  the  same  fixed  site, 
each  having  elevation  angle  k^  and  fix  and  data  points  for  satellite  pass 
No.  i are  plotted  in  a rectangular  coordinate  system  at  points  Y^  = f(kj), 

Xj  = L^,  then  the  points  Y^,  X^  ideally  fall  on  a straight  line  for  which  the 
slope  is  ^ = ~ ; or  the  inverse  slope  giving  the  antenna  height  error  and 
the  value  for  L at  y = 0 will  give  the  correct  longitude  position.  This  is 
essentially  the  approach  to  obtain  a three-dimensional  fix  (Anderle,  1971; 
Greely,  1971).  As  can  be  seen  from  Figure  1,  the  direction  from  the  true 
fix  (L)  toward  the  SAT  NAV  computed  longitude  (L^)  is  the  opposite  from  the 
direction  of  the  passing  satellite  if  the  initialized  antenna  height  is  too 
high  (that  is  toward  the  W (position  Li)  for  an  E-passing  satellite),  or 
toward  the  satellite  if  the  initialized  height  is  too  low. 
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However,  the  writer  previously  found  (Berg,  1975,  and  in  press)  that 
with  operational  receivers  under  field  conditions  in  Colombia  and  on  board 
Kana  Keoki , better  results  (less  spread  in  the  longitude  position  fix)  are 
obtained  when  the  height  correction  and  longitude  fix  are  taken  only  for 
individual  "double  passes",  that  is,  for  the  same  satellite  at  two  consecu- 
tive orbits,  one  east  and  the  other  west  of  the  site.  For  the  final  fix  and 
height  the  mean  of  all  such  double  passes  for  longitude  and  height  are  taken. 

As  an  example,  take  the  first  satellite  double  pass  (No.  67)  in  Table  1A,  B 
HNL  site  No.  25,  that  is  pass  3 and  4 (on  day  307).  (Refer  to  Fig.  2 for  the 
following  discussion.)  The  first  one  (No.  3)  was  E of  the  site  at  a 
K = 18°  (f(18°E)  = 0.000381),  and  the  second  (No.  4)  was  W of  the  site  at  a 
K = 31°  (f (31°W)  = - 0.000512).  The  longitude  of  the  first  one  was  -157°52.023' 
and  that  of  the  second  was  -157°52 .033 ' . The  latitude  for  both  (to  sufficient 
accuracy  for  the  height  correction)  was  21°18.8'.  The  antenna  height 
correction  is  then  given  by 

dHfml  = dx(m*-n  arc)  • COS(IAT) 
m dy (nautical  mile/m) 

(since  nautical  miles  = min  arc  • COS(IAT));  the  numerical  values  thus  give 


dH 


L1~L2 

f(k1)-f(k2) 


_ "52.023 -(-52. 033) 

cos  (LAT)-  0.000381-(-0. 000512) 


cos(21°18.8') 


= + 10.43  m. 
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Adding  the  correction  dH  = 10.43  m to  the  initialized  antenna  height  (23  m) 
then  yields  the  ANT  HT  of  33.43  m listed  in  Table  IB  (pair  1).  Since  the 
initialized  antenna  height  was  too  low  the  true  longitude  is  W of  the  first 
fix  (Lp  and  the  fix  error  is  dH  • f (kj) /cos (LAT)  = 0.0042',  so  that  the 

true  fix  for  this  double  pass  is  157°52.023'  + 0.0042'  = 157°52 .0272 'W  (also 

listed  in  Table  IB  (pair  1)). 

In  the  actual  computer  output  plots,  in  the  Appendix  as  well  as  in 
Figure  2,  the  increasing  absolute  value  of  the  longitude  is  always  in  the 
positive  x direction;  E-passing  satellites  are  plotted  positive  in  y(f(k)) 
if  longitude  is  E and  are  negative  in  y if  longitude  is  W.  As  a result,  the 
slope  of  the  line  connecting  the  two  passes  in  the  x,  y diagram  is  always 

positive  for  a positive  (upward)  correction  of  antenna  height.  The  more 

near  vertical  the  connecting  lines  the  less  antenna  height  correction  is 
indicated.  Note  that  the  x scale  is  different  from  plot  to  plot,  depending 
or.  the  range  in  longitude  of  the  single  pass  fixes  actually  observed  at  the 
site  and  used  for  the  double  pass  determinations.  On  each  computer  output 
plot  of  the  appendix,  the  two  positions  belonging  to  a double  pass  are 
connected  by  a line.  Above  the  line  is  a number  indicating  the  satellite 
identification  followed  by  a letter  N or  S (if  available)  indicating  that 
the  direction  of  satellite  travel  is  north-  or  south-bound.  Below  the 
connecting  line  the  number  indicates  the  chronological  order  of  the  double 
pass  in  a set  (that  is  the  pair  number  in  the  example  of  Table  IB  an.  :he 
corresponding  output  plot  in  the  Appendix,  Group  I,  site  HNL  25).  All 
double  pass  height  and  longitude  determinations  are  given  in  separate 
computer  printouts  (similar  to  Table  1A  to  D)  where  the  mean  standard 
deviation  (o)  of  a single  double  pass  observation,  and  standard  deviation 
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of  the  mean  for  the  set  of  observations  (a  ) are  calculated.  If  rejection 

m 

criteria  have  been  applied  to  eliminate  some  observations  further  calculations 
are  executed  (see  later) . The  result  of  these  calculations  is  also  printed 
at  the  top  of  the  computer  output  plots  of  the  appendix.  The  bottom  line 
of  the  output  plots  gives  general  information  as  to  date,  location  and/or 
site  number,  and,  after  "HT:",  the  initialized  antenna  height  (in  meters). 

In  order  to  obtain  geoidal  height  from  antenna  height  the  latter's 
distance  above  water  has  to  be  subtracted.  On  the  Glomar  Challenger  this 
antenna  height  is  near  22  m and  on  the  Kana  Keoki  it  is  18  m. 

Selection  criteria 

As  the  programs  were  developed,  it  was  found  useful  to  limit  the  double 
pass  fixes  to  those  judged  most  reliable.  Two  types  of  criteria  have  been 
used:  the  first  type  selects  satellite  passes  on  the  basis  of  information 

supplied  by  satellite  navigator  output,  and  the  second  selects  those  double 
pass  results  most  consistent  among  themselves,  or  rejects  "bad"  ones. 

The  first  type  of  criteria  that  have  been  used  throughout  involves 
retaining  satellites  with  elevation  angles  between  5°  and  75°,  and  their 
fixes  if  the  navigator's  computer  calculations  show  five  interactions  or 
fewer.  Another  option  uses  only  those  fixes  obtained  with  a prescribed 
minimum  number  of  symmetric  Doppler  counts  on  each  side  of  the  point  of 
closest  approach  (for  HNL  site  25  see  Table  IB,  top).  Selection  of  double 
pass  fixes  for  a single  day  or  by  the  satellite  identification  number  is 
possible.  For  all  Kana  Keoki  data  the  same  identification  numbers  as  in 
Daugherty  (1974)  or  Berg  (1975)  have  been  used.  Satellites  used  for  the 
DSDP  site  fixes  are  identified  by  a two-digit  code  corresponding  to  the 
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third  and  fourth  digits  in  the  Navy  satellite  number  in  the  following 
manner: 


Navy  Satellite 
Number 

DSDP  ID 
Number 

Kana  Keoki 

30  120 

12 

42 

30  140 

14 

54 

30  180 

18 

63 

30  130 

13 

64 

30  190 

19 

65 

30  200 

20 

99 

The  rejection  criterion  relates  to  the  standard  deviation  (a)  of  the 
double  passes  in  a set  of  observations,  that  have  be  retained  after  the 
selection  criteria  had  been  applied.  Since  erratic  or  bad  fixes  are  obtained, 
it  is  desirable  to  eliminate  those  from  the  data  set.  Rejection  of  positions 
outside  a 3a  (or  2cr)  limit  works  well  on  relatively  large  data  sets  (>  30 
double  passes)  , whereas  even  a 2cr  limit  "misses"  where  just  a quick  look 
would  subjectively  reject  one  double  pass  in  relatively  small  data  sets 
(four  to  ten  double  passes).  As  an  example,  see  the  output  plot  of  DSDP 
Site  353  (Fig.  3),  where  clearly  the  position  at  the  "0"  crossing  (longitude) 
and  the  inclination  of  the  connection  line  for  satellite  14/5  and  13/4  are 
different  from  the  others,  but  were  not  eliminated  by  the  2a  limit  applied 
in  the  rejection  criterion.  For  the  purpose  of  this  report  the  rejection 
with  a prescribed  multiple-of-a  limit  on  both  sides  of  the  mean  was  first 
applied  to  the  double  pass  results  of  the  antenna  height.  If  one  or  more 
heights  were  found  outside  the  prescribed  limits,  the  double  passes  were 
completely  dropped  from  the  list,  means  and  standard  deviations  were 
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recalculated  and  the  same  process  was  iterated  until  no  values  remained 
outside  the  na  limit  of  the  last  step  (see  example  HNL  site  25,  2a  iteration 
in  Table  1C) . After  elimination  of  all  double  passes  by  the  iteration  on 
height,  the  remaining  double  passes  were  iterated  in  a similar  manner  but 
as  to  longitude  (example  Table  ID).  Usually  the  bad  values  disappear  with 
the  first  few  iterations  on  height.  After  the  iteration  on  height  and 
longitude  are  completed,  only  double  pass  fixes  that  have  survived  these 
rejections  remain  for  calculation  of  mean  height,  longitude  (and  latitude). 

For  the  example  of  HNL  site  25  the  final  mean  height,  mean  longitude,  and 
mean  latitude,  together  with  the  a and  am,  are  printed  at  the  bottom  of 
Table  ID.  The  same  values  appear  on  the  top  portion  of  the  corresponding 
computer  output  plot  in  the  Appendix  Group  1 with  the  bottom  line  "DAY  306-313 
1972  KK  HNL  HT:  +23M  #25  2 ct",  where  25  is  the  site  number  (see  Appendix, 

description).  From  a different  example,  Figure  A shows  the  distribution  of 
height  values  for  HNL  site  41A  (that  is  the  1973  portion  of  site  41)  remaining 
after  a 3s  iteration  rejection.  The  four  values  shown  with  open  circles  are 
those  that  have  been  rejected  on  the  basis  of  subsequent  3a  iterations  on  the 
longitudes.  Figure  5 shows  the  longitude  distribution  of  the  same  site, 
remaining  after  the  3a  height  iterations.  Values  shown  by  open  circles  again 
indicate  those  rejected  by  the  3a  iteration  on  longitude. 

One  might  consider  proceeding  with  some  kind  of  absolute  limit  on  the 
excursion  from  the  mean.  Knowing  that  the  accuracy  of  the  NAV  SAT  system  is 
of  the  order  of  a few  tens  of  meters,  one  could  eliminate,  say,  the  two  to 
five  values  the  farthest  from  the  mean  (for  a larger  data  set)  and  additionally 
reject  all  height  and  longitudes  outside  of,  say,  70  m of  the  mean  (after 
elimination  of  the  first  two  to  five  extreme  values).  It  was  found,  however, 


that  proceeding  on  a ns  basis  allows  better  recognition  of  the  quality  of  a 
data  set  (as  given  by  a,  the  standard  deviation  of  a double  pass  fix).  A 
reasonably  good  fix  is  obtained  with  a about  20  m or  less  in  each  height  and 
longitude.  Since  height  and  longitude  are  not  independently  determined,  a 
height  accuracy  in  one  usually  is  associated  with  a poorer  determination  in 

f~~2  r1 

the  other.  Typical  values  for  y + a (obtained  for  HNL  site  41)  are 
26.7  n for  3a  iterations  and  12.8  m for  2a  iterations,  where  of  the  130 
double  passes  in  the  list  117  remained  after  3a  iterations  and  65  remained 
after  2a  iterations.  Smaller  data  sets  show  somewhat  larger  values. 

To  avoid  confusion,  one  limitation  of  the  present  computer  program 
should  be  pointed  out.  If  satellite  fixes  were  taken  at  a site  but  the 
initialized  antenna  height  was  changed,  the  program  runs  as  a data  set  only 
those  fixes  for  which  the  initialized  antenna  height  was  identical.  There 
were  instances  of  initialized  antenna  height  changes  on  both  the  Kana  Keoki 
and  the  Glomar  Challenger. 


ACCURACY  AND  RELIABILITY 

Stansell  (1976,  p.  481)  has  indicated  that  the  root  sum  square  error  of  a 
position  fix  is  in  the  range  of  18  to  35  m when  the  laboratory  standards  and 
sophisticated  refraction  correction  models  of  the  Applied  Physics  Laboratory 
are  used,  whereas  field  results  would  lie  in  the  range  of  27  to  37  meters.  In 
these  measurements,  the  exact  elevation  of  the  antenna  is  usually  known  and 
initialized.  In  the  operational  shipboard  SAT  NAV,  however,  only  a standard 
troposphere  is  assumed  for  the  fix  calculations.  Additional  errors  that  are  pre- 
dominantly reflected  in  the  longitude  fix  result  from  the  error  in  the  initialized 
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antenna  height.  It  is  this  later  error  that  is  reduced  by  the  double  pass 
method;  at  the  same  time  the  method  yields  the  antenna  height,  a result  not 
given  by  the  on-board  satellite  navigator  (for  a fixed  position). 

It  was  somewhat  of  a surprise  that  the  shipboard  results,  using  double 

passes,  indicated  similar  accuracies  in  height  and  longitude  as  those 

mentioned  by  Stansell,  especially  if  a larger  number  of  double  passes  are 

obtained.  With  the  exception  of  problem  sites  (like  Callao,  Peru)  the 

double  pass  method  gives  a standard  deviation  (for  one  double  pass  observation) 

of  less  than  20  m both  in  height  and  longitude  on  a 2 a iteration  for  ten  or 

more  double  pass  observations  where  all  individual  fixes  have  been  used  from 

satellites  with  elevation  angles  between  5°  and  75°,  and  five  or  less 

iterations  of  the  SAT  NAV,  where  no  restrictions  have  been  put  on  the  minimum 

Doppler  counts  (see  Table  2 and  Appendix  Group  I for  heights  determined  on 

board  Kana  Keoki) . Some  "good"  sites  even  yield  standard  deviations  of 

between  7 and  10  m on  a 2a  iteration  (KK  sites  14,  15,  25,  and  41  in  Table  2). 

Since  height  and  longitude  are  coupled  in  a double  pass  determination, 

however,  very  good  results  for  one  usually  lead  to  poorer  results  for  the 

2 2 

aL  +aH  , however,  is  usually  between 

25  and  30  meters . 

Since  in  the  years  1971  to  1974  the  Kana  Keoki  made  numerous  stops  in 
certain  harbors,  the  actual  position  in  each  harbor  being  different  from 
occasion  to  occasion,  we  restricted  the  tables  (Tables  3 and  4*)  to  comparison 
of  the  height  determination.  It  should  be  mentioned  that  no  (height) 


*Whenever  there  appears  a reference  to  a table  the  corresponding  computer 
output  plots  are  found  in  the  appendix.  Refer  to  the  description  of  the 
appendix  for  locating  data. 


corrections  for  ocean  or  solid  earth  tidal  variations  have  been  applied 


Such  variations  should  not  affect  the  double  pass  determinations  since 


changes  are  rather  small  over  the  time  span  of  1 hour  40  min  to  1 hour 


50  min 


in  the  locations  considered  here,  but  could  affect  the  a where  the  tidal 


ranges  are  several  meters.  In  harbors  like  Anchorage,  Alaska  with  some  10  m 


changes  the  tidal  variations  should  be  taken  into  account 


bad  values  (in  height  and  longitude).  Bad  values  occur  sometimes  for 


individual  fixes  when  the  total  number  of  Doppler  counts  is  rather  small  and 


very  asymmetric,  the  iterations  required  by  the  SAT  NAV  are  four  or  five 


but  sometimes  we  were  unable  to  determine  any  reason  based  on  returned  data 


the  computer  does  not  always  catch  a bad  value  that  could  be  eliminated  by 


eye  or  perhaps  in  a more  sophisticated  manner  by  specifying  an  absolute  value 


of  the  deviation  from  the  mean 


Variation  with  time  of  da 


Ihe  reason  for  applying  2 ct  iterations  to  the  larger  data  sets  is  somewhat 


Some  more  extreme  values  occur  at  given  (local)  times  of  the  day 


probably  due  to  changes  in  the  tropospheric  refraction  index.  The  author 


previously  found  (Berg,  1975)  that  double  pass  longitudes  showed  less  spread 


in  the  late  night-early  morning  hours  for  Suva.  The  2a  iterations  tend  to 


eliminate  the  more  extreme  values.  Figure  6 shows  the  height  distribution 


for  one  of  the  largest  data  sets  (HNL  41)  as  a function  of  time  (in  GMT) 


0 hrs  GMT  is  14  hrs  local  in  HNL.  There  is  a well-marked  trend  of  a steadily 


decreasing  antenna  height  between  1 hrs  GMT  to  7 hrs  30  min  and  a rather 
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snall  variation  between  9000  and  1730.  The  times  as  given  here  for  a double 
pass  are  those  of  closest  approach  of  the  first  pass  of  the  satellite;  the 
second  pass  occurred  some  1 hr  46  min  later.  The  antenna  height  determined 
from  the  46  values  between  9000  and  1730  hr  is  19.1  m ± 8.4(a)  and 
- 1.24(aa),  which  comparesto  20.96  ± 8.94(a)  ± l.ll(am)  determined  by  the 
2a  iterations  on  site  41  (Table  4A) . The  value  of  19.1  m is  probably 
somewhat  low  because  of  the  three  low  points  in  Figure  6.  The  even  more 
time-restricted  set  between  1300  and  1518  yields  22.90  ± 5.37(a)  ± 1.39(am)  m 
from  only  15  double  pass  observations,  indicating  that  the  local  night  hours 
from  about  3 a.m.  to  near  7 a.m.  (time  of  the  last  second  pass)  give  a much 


more  reliable  observation  and  that  the  am  is  only  slightly  larger  than  that 
of  the  total  set.  The  night  observations  (22.9  m)  can  be  compared  to  what 


the  writer  thinks  is  the  best  height  determination  for  Honolulu,  that  is,  a 

2a  iteration  on  heights  alone  (hand -calculated)  which  retains  198  from  all 

248  heights  determined  from  1972  to  1974  (see  Table  4A)  of  22.58  ± 13.57(a) 

- 0.96(aa)  meters.  Perhaps  taking  the  ocean  tides  (about  .6  meter)  into 

account  would  further  reduce  the  standard  deviation  of  early  morning 

observations  to  less  than  5.0  meters. 

Figure  7 similarly  shows  the  longitude  double  pass  fix  distinctions 

as  a function  of  time  of  day  for  HNL  41A. 

Table  4A  (for  HNL)  and  Table  3A  (for  Suva)  further  compare  height 

results  at  different  port  calls  and/or  locations  in  the  harbor.  The  a 

m 

(standard  deviation  of  the  mean)  limits  usually  straddle  the  best  determina- 
tions that  have  been  obtained  by  2a  iterations  (on  the  height  alone)  from 
all  heights  obtained.  To  be  more  precise,  the  confidence  intervals  should 
be  determined  for  each  data  set.  As  a function  of  the  number  of  observations 
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the  95 7,  (98%)  confidence  intervals  in  multiples  of  om  are:  for  four 

observations  3.2am  (4.5a  ) and  for  120  observations  1.98a„  (2.36a  ) (see 
Crow  et  al ■ . 1960,  p.  231).  The  best  Honolulu  determination  of  22.6  m has 
a 95%  (987,)  confidence  interval  of  ± 1.9  (2.3)  m,  and  the  early  morning 
value  of  22.9  m has  a 957,  (987,)  confidence  interval  (15  observations)  of 
± 3.0  (3.7)  meters. 

A different  way  of  looking  at  this  accuracy  problem  consists  of 
determining  the  convergence  of  the  mean  values  with  increasing  number  of 
observations.  This  method  is  illustrated  in  Figure  8 for  the  height  and 
longitudes  that  remained  for  Honolulu  site  41A  after  the  3a  iterations 
processing. 

The  first  data  point  corresponds  to  the  first  double  pass  height  (and 
longitude)  determination.  The  next  point  (N  = 2)  is  the  mean  value  obtained 
for  the  height  (longitude)  from  double  pass  one  and  double  pass  two.  The 
N'th  point  is  the  mean  after  N double  passes.  Nearly  stationary  values  are 
obtained  after  30  double  passes.  Oscillations  in  antenna  height  are  clearly 
associated  with  time  of  day,  the  higher  peaks  being  associated  with  the 
early  hours  (in  GMT)  of  each  day,  as  they  are  in  Figure  6.  In  addition,  the 
final  values  show  the  ± lam  interval  for  this  3 a iteration  and  compare  them 
to  (a)  the  final  value  obtained  through  2a  iterations;  (b)  the  final  values 
obtained  for  the  full  set  for  site  41  for  both  the  3a  and  the  2a  iterations 
(the  other  selection  criteria  are  identical)  and  (c)  the  heights  of  the  3a 
and  2a  Iteration  values  on  all  heights  obtained  for  HNL.  As  can  be  seen 
from  Figure  8,  the  2 a iterations  on  the  different  sets  give  somewhat  lower 
heights  than  the  3a  iterations,  since  the  former  tend  to  eliminate  the 
consistently  higher  values  observed  in  the  early  day  time  (GMT)  hours.  In 


Daugherty  (1974)  has  mentioned  that  the  reason  for  the  small  standard 
deviation  on  longitude  and  latitude  for  Wellington  (site  15  in  Table  2)  was 
the  high  number  of  symmetric  Doppler  counts.  We  have  therefore  analyzed  in 
some  detail  the  data  for  a number  of  sites  to  compare  the  dependence  of  the 
double  pass  results  on  increasing  number  of  symmetric  Doppler  counts.  In 
the  tables  and  output  graphs  (Appendix  Group  III)  the  "number  of  symmetric 
Doppler  Counts"  or  "D.C.  > No"  means  that  there  have  been  more  counts  both 
before  and  after  the  closest  approach  than  the  number  indicated.  Thus  in 
Table  43  for  HNL  sites  40  and  41  in  the  last  line  under  the  heading  "SYM.CTS." 
the  number  ">  10"  means  that  there  have  been  at  least  11  Doppler  counts  before 
closest  approach  and  at  least  11  Doppler  counts  after  closest  approach. 

In  Table  4B  for  HNL  sites  40  and  41  the  antenna  heights  obtained  by  the 
3a  and  2a  iterations  on  all  fixes  are  compared  to  those  obtained  by  restricting 
the  observations  to  progressively  higher  numbers  of  symmetric  counts  (and 
using  the  liberal  3a  iterations  to  eliminate  bad  fixes).  As  can  be  seen  in 
Table  4B  and  Appendix  Group  III  the  height  determination  becomes  better  (the 
standard  deviation  diminishes)  for  increasing  number  of  symmetric  counts. 

But  even  for  fixes  obtained  by  using  only  passes  with  a high  number  of 
symmetric  Doppler  counts,  there  remain  some  (in  the  3a  iterations)  that  one 
would  like  to  eliminate.  Where  they  actually  were  eliminated,  the  results 
are  indicated  by  an  asterisk  in  Table  4B. 
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Tables  3B  and  4B  indicate  that  little  final  accuracy  (cj^)  is  gained  by 

restricting  double  pass  height  determinations  to  those  with  a high  number  of 

symmetric  Doppler  counts,  that  the  final  height  values  are  close  to  those  of 

the  2a  iterations  on  all  data,  and  that  because  of  the  larger  amount  of  fixes 

available,  the  2a  iterations  have  smaller  standard  deviations  of  the  mean 

(a  ) and  smaller  confidence  intervals, 
m 

It  should  be  pointed  out  that  even  highly  symmetric  Doppler  counts 
sometimes  lead  to  fixes  that  are  off  from  the  remaining  data  as  was  the  case 
for  one  (single  fix)  with  27  total  Doppler  counts  and  13  after  closest 
approach;  elevation  angle  was  22°  and  the  satellite  navigator  did  only  two 
iterations  for  the  fix.  This  questionable  fix  leads  to  an  antenna  height  of 
118  m for  the  double  pass  determination  in  the  HNL  40  set  and  compares  to  the 
12.7  m for  the  remaining  seven  double  passes  (in  the  last  line  of  Table  4B) 
for  symmetric  Doppler  counts  > 10.  Eliminating  this  bad  height  also  reduced 
the  a from  40.3  to  16.7  m,  more  in  line  with  the  accuracy  for  double  passes. 
That  the  height  value  (12.7  m)  of  this  particular  set  of  double  passes  is 
lower  than  the  usual  22  m for  HNL  stems  from  the  bias  due  to  a particular 
satellite  (65)  that  consistently  seems  to  indicate  lower  heights  (Berg,  1975) . 
Results  for  additional  sites  obtained  by  selecting  higher  numbers  of  symmetric 
Doppler  counts  are  presented  in  Table  5 and  are  more  extensively  discussed 
under  "Problem  Sites". 

Variations  among  different  satellites 

In  the  earlier  report  (Berg,  1975),  no  significant  differences  were  found 
I 


X 


between  the  double  pass  fixes  of  different  satellites  for  the  longitudes.  It 
was  noted,  however,  (Berg,  1975)  that  satellite  65  always  showed  a lower 
antenna  height  than  the  other  ones.  Si.ee  in  the  earlier  report  longitudes 


> 


-19- 

and  heights  were  obtained  graphically  and  from  somewhat  different  numerical 
values  for  the  height  corrections,  the  Suva  data  (sites  13  and  14)  have  been 


» 


t 


» 


I 


reexamined  together  with  those  for  HNL  (site  41)  where  a sufficient  number  of 
double  passes  for  individual  satellites  were  available. 

As  can  be  seen  from  Table  3C  for  Suva  (sites  13  and  14)  and  Table  4C  for 
HNL  (site  41)  (Appendix  Group  IV),  the  a of  a height  determination  from  a 
specified  satellite  is  of  the  same  order  or  better  as  that  from  all  satellites 
combined . 

Combining  all  heights  determined  for  satellites  42,  54,  63  and  64  and  65, 
however,  a value  is  obtained  (59.5  m)  different  from  that  of  65  alone  (35.8  m) 
for  Suva  sites  13  and  14.  Similarly  the  data  for  HNL  site  41  indicate  (Table  4C 
and  Fig.  9)  a lower  height  (10.5  m for  2a  iterations)  for  satellite  65  than  for 
all  satellites  combined  (21.0  m for  2a  iterations).  This  lower  height  also  is 
obtained  when  only  symmetric  Doppler  counts  > 10  are  selected  (together  with  3a 
iterations).  There  are  four  such  double  pass  determinations  for  HNL  site  40  with 
a height  of  1.57  ± 5.52(a)  ± 2.76(am)  m and  six  such  double  pass  determinations 
for  HNL  site  41  with  a height  of  6.50  ± 15.64(a)  ± 6.39(am)  meters. 

The  fact  that  the  height  difference  between  satellite  65  and  the  remaining 
ones  is  larger  in  Suva  (in  1971)  than  in  Honolulu  (in  1973/1974)  may  indicate 
that  the  height  is  either  site-  or  epoch-dependent.  Perhaps  the  accuracy  of 
the  orbital  parameters  has  increased  with  time  through  better  modelling  of 
the  earth's  gravity  field. 

Problem  sites 

Problem  sites  are  those  where  the  fixes  returned  by  the  satellite 
navigator  show  a rather  large  scatter  when  the  ship  is  in  a fixed  position, 
a sufficient  number  of  satellite  passes  is  available,  and  the  initialized 
antenna  height  is  sufficiently  close  to  the  true  height.  Callao  (Peru)  was 
recognized  early,  even  from  on  board  ship,  as  being  such  a site.  Using  the 


J 


-20- 

double-pass  method  and  3a  or  2ct  iterations  on  the  data  does  not  eliminate 
the  large  spread,  and  a problem  site  therefore  can  be  considered  more 
specifically  as  one  where  the  a of  the  double-pass  observation  is  larger 
than  on  other  sites,  especially  on  the  antenna  heights  for  a sufficiently 
large  data  set.  Table  1 and  the  corresponding  computer  plots  (Appendix 
Group  I)  show  that  even  after  the  tight  2a  iterations,  the  standard  deviation 
of  single  double-pass  height  for  Callao  in  1972  and  1974  (as  well  as  for 
Acapulco  in  1972)  gives  rather  large  values  or  more  than  40  meters.  A 
similar  large  standard  deviation  is  obtained  for  the  double-pass  longitude 
determinations  (see  computer  plots.  Appendix  Group  I,  for  sites  43,  44,  and 
45).  The  data  in  1972  seem,  however,  of  better  quality  (sites  16  and  17). 

As  Table  5 indicates,  selection  of  double  passes  with  a high  number  of 
symmetric  Doppler  counts  (>  8)  does  not  reduce  the  spread.  The  particularly 
bad  data  on  site  44  (Callao  B)  might  be  due  in  part  to  the  fact  that  a 
large  grain  elevator  was  quite  close  to  the  docking  position  of  the 
Kana  Keoki  so  the  selection  by  symmetric  Doppler  counts  was  not  run  on  the 
computer.  Similarly  the  longitudes  of  the  double-pass  fixes  for  site  43 
(15  observations  remaining  with  > 8 symmetric  counts  and  3a  iterations) 

I 

still  show  a large  spread  and  a large  76-m  standard  deviation  (see  computer 
plots) . 

The  writer  was  told  by  Hawaii  Institute  of  Geophysics  personnel  that 
once  the  ship  leaves  port,  the  fixes  obtained  on  route  seem  normal  again. 

The  reasons  for  the  larger  spread  of  fixes  in  Callao  are  not  known.  Radio 
(or  television)  interference  might  cause  problems,  or  atmospheric  conditions 
might  cause  anomalously  large  changes  in  the  tropospheric  refractive 


index. 


SEA  LEVEL  AND  GEOIDAL  HEIGHTS 


As  stated  earlier,  reasonably  good  results  are  obtained  when  the 
standard  deviation  (ct)  of  the  antenna  height  or  longitude  of  a single 
double -pass  observation  in  a larger  data  set  is  20  m or  less  after  the  2 ct 
iterations.  The  fixed  position  requirement  is  certainly  fulfilled  for  the 
fixes  obtained  by  the  R.  V.  Kana  Keoki  when  it  is  tied  to  a dock  in  a harbor 
and  good  results  were  also  obtained  in  one  case  when  the  ship  was  anchored 
(see  site  21,  Punta  Arenas,  Table  2 and  the  corresponding  computer  plot 
Group  1) . 

As  for  the  Glomar  Challenger,  the  position  of  the  ship  between  two 
consecutive  orbits  of  a double  pass  does  not  seem  to  vary  enough  to  affect 
the  antenna  height  results.  In  the  case  of  DSDP  site  410/4 10A,  the  longitudes 
of  the  SAT  NAV  output  were  corrected  to  the  hydrophone  beacon  position 
(assuming  x = 0,  y = 0 is  at  the  location  of  the  drill).  The  ship's  bow 
direction  (from  the  gyro  heading)  was  also  taken  into  account.  These 
corrections  were  applied  to  double  passes  of  satellite  No.  12.  The  antenna 
height  found  for  the  five  double  passes  of  satellite  12  (retained  after  the 
2ct  iterations  of  the  computer)  was  74.13  m (ct  = 5.48)  as  compared  to  74.89  m 
(ct  = 5.21)  when  position  corrections  were  applied  and  75.79  m (ct  = 7.30) 
from  all  15  double  passes  retained,  but  uncorrected. 

A better  (satellite-determined)  beacon  position  is  obtained  for 
longitude  and  latitude.  The  ct  = 25.8  m for  the  uncorrected  longitude  values 
of  the  five  double  passes  retained  of  satellite  12  reduces  to  18.4  m after 
corrections  and  compares  to  27.1  m for  all  15  uncorrected  double  passes. 
Similarly  applying  corrections  to  the  latitude  data  (satellite  12,  five 
double  passes)  reduces  the  cr  by  107».  It  would  appear  therefore  that  first 
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reducing  the  SAT  NAV  fixes  to  the  beacon  position  would  improve  the  overall 
accuracy  of  the  site  coordinates.  The  Glomar  Challenger  double  pass  fixes 
at  drill  site  location  seem  at  least  as  accurate  as,  if  not  better  than,  the 
harbor  data  on  board  the  Kana  Keoki . This  might  be  due  to  the  difference  in 
epoch,  to  the  DSDP  data  being  more  recent,  to  the  incidentally  better 
elevation  angle  distribution  resulting  in  a higher  number  of  symmetric 
Doppler  counts  and  to  more  stable  tropospheric  conditions.  In  retrospect 
it  also  appears  that  the  very  slow  speeds  involved  in  moving  the  Glomar 
Challenger  around  the  beacon  position  (of  the  order  of  a few  0.001  knots), 
do  not  influence  the  fix  results.  This  conclusion  seems  reasonable  when  one 
compares  the  fix  errors  resulting  from  speed  errors  in  northerly  or  easterly 
direction  (Stansell,  1976,  his  figures  p.  490  and  491)  over  most  of  the 
satellite  maximum  elevation  angle  range  used  for  the  double  pass  determinations. 

Finally  in  Table  6 and  Figure  10,  sea  level  height  with  respect  to  the 
Navy  navigation  satellite  reference  ellipsoid  are  compared  to  two  different 
geoid  models.  The  sea  level  heights  are  those  determined  from  double  passes 
during  the  present  investigation  by  subtracting  the  antenna's  height  above 
the  ship's  waterline  from  the  antenna  height  determined  from  the  double 
passes.  Details  of  antenna  height  have  been  supplied  by  Ted  Gustafson 
(Scripps  Institution  of  Oceanography)  to  the  nearest  inch  for  DSDP  sites  407, 
410,  and  412.  For  other  DSDP  sites  this  height  was  assumed  to  be  21.9  m 
(last  column  in  Table  6).  The  height  for  the  Kana  Keoki  antenna  was  taken 
from  the  ship's  plans  and  a waterline  when  the  ship  was  in  Honolulu  and  is  close 
to  18  meters.  Geoidal  heights  were  read  off  two  different  maps.  The  one 
used  under  "SAT  NAV  MANUAL"  is  identical  to  the  one  published  by  Stansell 
(1976,  his  figure  p.  488)  with  10-m  contour  intervals  (missing  in  some  parts 
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of  the  map) . The  map  used  under  GEM  6 was  a model  prepared  and  published 
by  the  Defense  Mapping  Agency  Topographic  Center,  Washington,  D.  C., 

June  1974.  Figure  11  shows  the  Kana  Keoki  sites  and  DSDP  site  together 
with  GEM  6 geoidal  heights.  Since  both  maps  were  rather  small,  geoidal 
heights  in  Table  6 are  probably  no  better  than  1 to  2 meters  near  smooth 
topography  and  errors  are  possibly  larger  at  locations  with  steep  gradients. 
The  standard  deviation  (cr)  and  the  deviation  of  the  mean  (am)  in  Table  6 
are  those  of  the  data  remaining  after  the  (usually  2a)  iterations  on  height, 
followed  by  iterations  on  longitude  of  the  remaining  data,  a and  am  are 
therefore  too  small.  As  a consequence  the  95%  confidence  interval  of  the 
mean  based  on  crm  and  the  number  of  remaining  observations  is  also  somewhat 
too  small  (Table  6,  Fig.  10). 

As  Table  6 and  Figure  10  show,  despite  the  underestimate  of  the  95% 
confidence  intervals,  the  sea  level  heights  determined  at  all  the  DSDP  sites 
in  the  open  ocean  are  within  the  95%  confidence  interval  when  compared  to 
the  GEM  6 geoidal  height  model  and  when  compared  to  the  geoidal  height  map 
of  the  SAT  NAV  manual,  except  for  site  354  where  a rather  large  difference 
exists  between  the  SAT  NAV  manual  and  GEM  6 model.  This  site  is  located 
north  of  the  north  coast  of  Brazil  in  the  geoidal  valley  leading  to  the 
minimum  near  60°W  longitude  and  20°N  latitude  on  the  GEM  6 map. 

Sea  level  heights  determined  in  harbors  onboard  the  Kana  Keoki  in  some 
locations  seem  to  differ  from  the  GEM  6 geoidal  height  model.  Among  those 
locations  are  (see  Table  6,  Fig.  10):  Suva  (data  taken  1971,  1972), 

Wellington  (1971),  Quayaquil  (1972)  and  Panama  (1974),  all  located  in  areas 
of  relatively  steep  geoidal  height  gradients. 


The  largest  deviations  from  the  GEM  6 model  at  Panama  and  Quayaquil 
are  near  20  meters.  As  Marsh  nt  al.  (1976)  have  shown,  differences  of  more 
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than  20  m exist  between  the  GEM  6 and  SAO  SE  III  models  and  also  between 
the  Sky  Lab  altimeter  data  from  a continuous  revolution  on  January  31,  1974 
and  th'*  GEM  6 model. 

The  GEM  8 model  map  published  by  King-Hele  (1976)  is  too  snail  and  the 
print  too  weak  to  allow  determination  of  the.  heights  of  a number  of  sites; 
however,  most  of  the  sea  levels  given  here  seem  to  agree  better  with  GEM  8 
than  GEM  6,  where  a reasonable  check  was  possible. 

CONCLUSIONS 

The  double  pass  method  allows  three-dimensional  fixes  from  the  Navy 
navigation  satellite  system  for  a fixed  antenna  position;  only  latitude, 
longitude,  satellite  elevation  angle,  geometry  of  the  passes,  and  the 
satellite  identification  are  needed  from  the  satellite  navigator  printout. 
Sea  level  or  geoidal  heights  therefore  are  obtained  from  shipboard  fixes  in 
harbors  and  from  the  Glomar  Challenger  when  it  is  in  a fixed  position  over 
an  acoustic  transponder  in  the  open  sea. 

The  accuracy,  measured  by  the  standard  deviation  (a)  of  a double  pass, 
is  20  m or  better  in  height  and  longitude  after  "erratic"  or  "bad"  data  have 
been  eliminated  that  are  outside  a specified  limit,  for  this  investigation 
usually  2rj  in  height,  followed  by  elimination  of  data  outside  2a  in 
longitude,  for  sufficiently  large  data  sets. 

The  small  and  slow  motions  of  the  Glomar  Challenger  over  the  trans- 
ponder do  not  seem  to  affect  significantly  the  height  determinations, 
whereas  correcting  longitude  and  latitudes  to  the  transponder  coordinates 
seems  to  diminish  the  standard  deviation  of  these  observations. 
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Sea  level  determinations  at  eight  DSDP  open  ocean  sites  and  nine  Kana  Keoki 
harbor  sites  (Table  6 and  Pig.  10)  are  within  about  20  m of  the  model 
GEM  6 geoidal  heights.  Most  of  the  observed  heights  seem  to  be  closer  to 
those  of  the  GEM  8 model;  however,  the  smallness  and  weak  print  of  the 
published  (King-Hele,  1976)  map  did  not  allow  a tabulation.  The  standard 
deviations  of  the  mean  height  for  a four-day  observation  period  are  less 
than  4 m,  reaching  as  low  as  1.9  meters.  This  limit  can  be  reduced  further 
by  obtaining  Doppler  fixes  only  from  satellites  in  a double  pass  configuration. 

Most  important  perhaps  is  the  fact  that  open-ocean  sea  level  is 
determined  by  an  independent  method,  and  when  compared  to  the  Sky  Lab 
altimetry  data,  is  not  affected  by  uncertainties  in  a single  orbit  position 
determination  but  rather  represents  an  average  over  many  passes  of  several 
satellites . 

Incidentally  better  positions  for  the  Glomar  Challenger  ((x,y)  = (0.0) 
beacon  position)  are  obtained  when  beacon  position  corrections  are  applied 
to  the  satellite  navigator  positions.  Using  satellite  12  (30120)  for 
site  410,  a reduction  of  a from  25  to  18  m in  longitude  and  about  10%  in 
latitude  was  found.  It  appears  therefore  that  after  application  of  the 
reductions,  the  open-sea  fixes  are  of  the  same  quality  as  those  of  the 
Kana  Keoki  in  a fixed  harbor  location,  indicating  that  the  (x,y)  = (0.0) 
position  error  (cm)  after  four  days  of  satellite  navigation  fixes  should  be 
no  more  than  10  m in  longitude  and  in  latitude. 


Institution  of  Oceanography  cooperated  by  providing  the  satellite  navigation 
data  (through  Dr.  G.  P.  Woollard) ; special  thanks  are  due  to  Dr.  Ted  Gustafson 
who  provided  the  details  of  the  antenna  height  above  the  load  water  line  for 


the  Gloraar  Challenger  as  well  as  the  description  of  the  beacon  positions  of 
^ the  ship. 

Special  recognition  is  due  to  Duncan  M.  Chesley,  who  developed  the 
computer  programs  and  plotting  routines  used  for  presentation  of  the  results. 
I Finally  the  author  expresses  deep  appreciation  to  Dr.  G.  P.  Woollard 

for  his  continuous  interest  and  critical  discussions. 
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1 57 

52. 0 1 8 W 
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S-W 

54 

3 12 

1 7 
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1 8. 801 N 

l 57 

5 l . 993* 

2 

19 
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55 

312 
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64 

38 

21 

1 f . 825N 

157 

52.  022* 

2 

24 

7 
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56 

312 
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4? 

55 
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1 8. 835N 

157 

52. 025 a 

O 

A- 
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157 

5 ? . 0 1 6 a 
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S-W 

58 

312 

21 
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1 57 
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1 57 
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1 57 
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13 

0 
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2 
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6 3 

57 
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1 ri. 772N 
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5 2.  0 20  a 

2 
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c 

N-  A 
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SELECTION  CRITERIA 


PASSES  WITH  THt  FOLLOWING  TRAITS  HAVE  BEEN  RETAINED 


M 


»• 

FLEVATICN  ANGLE: 

5 TU  7 5 

DEGREES 

I TER  AT  I ONS: 

FEWER  T FAN  6 

SECOND  SIDE 

DOPPLER 

COLNT  : 

MORE  THAN 

I NDEX 

Sat 

PAIRS 

DAY 

1 

67 

3 

4 

30  7 

21 

2 

54 

1 3 

14 

307 

:* 

3 

63 

1 b 

17 

3 C 8 

4 

54 

19 

21 

308 

5 

65 

20 

22 

308 

22 

6 

63 

23 

24 

3 08 

:j 

7 

54 

27 

29 

308 

t* 

8 

65 

28 

30 

308 

9 

6 3 

31 

32 

309 

3’ 

10 

42 

33 

34 

309 

22 

1 I 

54 

36 

38 

309 

ii 

12 

65 

37 

39 

3 C 9 

13 

63 

41 

42 

3 C 9 

ii 

14 

42 

44 

46 

312 

Ja 

___  15  _ 

54 

48 

50 

312 

32 

16 

63 

51 

53 

312 

28 

1 7 

42 

54 

56 

312 

;* 

18 

64 

55 

57 

312 

•*<; 

1 0 

54  

58 

60 

312 

41 

20 

65 

59 

62 

312 

4? 

21 

63 

61 

63 

313 

I ) 


i 

1972  K<  HNL  HT 

: +23  M 

U 25 

6 

7 

i 

DATA  BELOW  ARE 

FOR  DAYS 

3C6-  313 

9 

U 

• l 

1 



. _ 

LONG 1TUDE 

4 

PAIR 

ANT.  HT. 

DEG. 

MINUTES 

*s 

I 

33.43 

157. 

52.0272 

7 

2 

37.34 

157. 

52.0341 

1 

3 

76.04 

157. 

51  .9765 

l 

4 

37.78 

157. 

52.0 1 77 

•J 

5 

4 1.36 

157. 

52.C229 

II 

6 

25.61 

157. 

52.0501 

. ? 

7 

26.71 

157. 

52.0194 

1 

8 

49. 70 

1 57. 

52 .0 1 93 

4 

9 

54.54 

157. 

52. C 4 46 

10 

36.15 

157. 

51  .8853 

1 1 

1 8.  29 

157. 

52.0344 

12 

26.23 

157. 

52.0  2 2 1 

1 

13 

15.57 

157. 

52.0 1 53 

.■» 

14 

42.62 

157. 

52.0129 

'> 

15 

27.63 

157. 

52.0342 

11 

16 

24.49 

157. 

52.0166 

1 7 

4 1.86 

157. 

5 2 • C C 4 8 

s* 

ia 

19.40 

157. 

52.0196 

'% 

1 9 

2 1.05 

157. 

52.0260 

■ » 

20 

26.31 

157. 

52.0219 

i 

21 

39  .20 

157. 

52.0016 

TABL2  1 B 


157?  KK  HNL  HT:  ♦ 23  M *25 


TABL 


J AT  A PEL  Li*  AS  fc  FUR  U AYS  30o-  313 


1 c 


MEAN  LCNGI T JOE  : 


15  7. 


52 . 3 1 A 5 MINUTES 


STD. 
DEV.  CE- 


DE V • 1 
mean: 


0.C334  MINUTES  l 57.70  METERS) 

C.0C73  MINUTES  ( 12.59  METERS > 


MEAN  latitude: 


21  . 


1 8 . 3 C A 2 MINUTES 


STD.  DEV.: 
DEV.  OF  mean: 


0.0336  MINUTES  ( 71. A6  METERS) 

0.006C  MINUTES  ( 11.03  METERS) 


MEAN  ANTENNA  HEIGHT; 
STD.  DEV.: 

OEV.  CF  mean: 


3A.3  5 METERS 
14. ie  METERS 
3.09  METERS 


***  LATITUDE  CORRECTION  APPLIED  *** 


. ....  ....  CORRECTING  HEIGHTS 

throut:  the  following  elements  have  been  removed: 

ELEMENT  3 WAS  2.9A  SIGMA  FROM  THE  MEAN 


MEAN  LONGITUDE: 

157. 

£2.01 66 

MI  NUTES 

STD. 

DEV  .: 

0 .0331 

MINUTES  I 

57.1  3 

ME  TERS  ) 

DEV . OF 

mean: 

0.C074 

MINUTES  ( 

12.78 

ME  TERS  » 

MEAN  LATI TUDE : 

2 1 . 

18.8054 

MINUTES 

STD. 

oev  .: 

0.0391 

MINUTES  ( 

72.55 

ME  TERS  ) 

DEV.  CF 

mean: 

C .0062 

MINUTES  ( 

11.47 

METERS ) 

MEAN  ANTENNA  HEIGHT:  32.26  METERS 
STD.  DEV.:  1C. 75  METERS 
DEV.  OF  mean:  2.  AO  METERS 


throut:  the  following  elements  have  been  removed: 

ELEMENT  9 WAS  2.07  SIGMA  FROM  THE  MEAN 


MEAN  LONG! TUDE : 


157. 


52.0151  MINUTES 


STD.  DEV.: 
DEV.  OF  mean: 


0.0  333  MINUTES  ( 
0.CC76  MINUTES  ( 


57.52  METERS) 
13.20  METERS) 


MEAN  LATITUDE: 


2 1 


13.8062  MINUTES 


STD.  DEV.: 
DEV.  OF  mean: 


0.0399 
C . 0 C 6 5 


MI  NUTc'S 
M I NUTES 


73.96  METERS) 
12.0  0 METERS) 


mean  antenna  height: 

STD.  DEV.: 

DEV.  OF  mean: 


31  . C9  ME  TERS 
9.65  METERS 
2.2  1 METERS 


throut:  no  values  outside  yuuk  2.00  sigma  limit 


heights  completed 


2 DROPPED  FROM  ORIGINAL  21 


• • . • . COMt'tL  T Ir.G  LUNUITUULS  • • • 


TA^LE  1 D 


THRCUTI  THf  f CiLLUa  INC*  EtLMt'  NTS  HAVE  Of:  f.N  HI  HOVLOS 
ELEMENT  10  -AS  J.VO  5 1 CM  A F ROM  iHt  MEAN 


MEAN  LONG! TUJC I 


STD.  utv.: 
DEV.  OF  mean: 


5?#  022  3 MINUTES 


0 .0  1 i 3 Ml  NOILS  ( 
C.CO  2/  MINUTES  ( 


19.67  MfTE«SI 
4.6  1 METERS! 


MF AN  L AT  I TUDE : 


sto.  ocv.: 
oev.  of  mean: 


IB.* I 23  MINUTES 


0 • C 223  Ml  notes  f 
C.C037  MINUTES  C 


41  .2  5 Ml  TEHSI 
6.BB  METERS! 


mean  antenna  height: 
stj.  oev.: 

DEV.  OF  mean: 


30.6  1 METERS 
9.85  VETERS 
2.32  Me.  TERS 


tmrcut:  t« e following  elements  have  oeen  removed: 


element  6 «AS  2.45  SIGMA  FROM  ThE  MEAN 


r _ MEAN  LONG1 TUOS: 


mean  l at i ruos: 


STC.  oev.: 
oev • or  mean: 


STO.  LEV.: 
DEV.  of  mean: 


52.02C7  MINUTES 


0.CO92  MINUTES  C 

0.0C22  minutes  ( 


16.6125  MINUTES 


0.0229  MINUTES  C 
0.CC39  MINUTES  ( 


MEAN  ANTENNA  HEIGHT: 
STO. • uEV. : 

DeV.  OF  Mg A HI 


31.1  2 METERS 
1C  .C  6 METERS 
2.4  4 METERS 


15. Q5  ME  TERS  1 
3.67  MF  TERS  f 


42.38  METERS) 
7.27  METERS) 


THPOUT:_THE  follc* ing  CLEMENTS  have  been  removed: 


ELEMENT  21  to  A 5 


2.07  SIGMA  FROM  THE  MEAN 


MEAN  LONG I TUDS  J 


STC.  dev.: 
DEV.  UF  mean: 


52.C2I9  MINUTES 


0.CC3I  MINUTES  C 
C.0C2C  MINUTES  ( 


13.95  METERS) 
3.49  ME  T £ RS  ) 


MEAN  LAf  XTUOc:  


STD.  C»:V.: 
dev.  of  mean: 


16.6142  MINUTES 


0.0222  MINUTES  ( 
C • 3 0 39  MINUTES  ( 


41.19  METERS) 
7.28  FILTERS! 


MEAN  ANTENNA  HEIGHT: 
STC.  CfV.: 

OE  '4  • Vr  MEAN: 


20.61  METERS 
1C  • 1 6 METERS 
2.54  METERS 


tnrout:  the  roll  )• ing  elements  have  been  remcveo: 
ELEMENT  17  0AS  2.11  SIGMA  FF CM  THE  MEAN 


me a n longitude: 


STD.  DEV.: 
DEV*  UF  mean; 


52.0230  MINUTES 


C.CC69  MINUTES  ( 
C *0C IB  M i MUTES  ( 


FINAL 


11.94  METERS!  RT?QTTTrTq 

3. CD  MCTCRS!  Hli3UU.3 


« M£AN  LA T I TUOE: _ . 


STD.  DtV.: 
OEV.  UF  mean: 


16.6140  MINUTES 


C.0275  MINUTES  ( 
0.C041  MINUTLS  I 


41.73  MfTCRS! 
7.62  MLTLRS) 


MEAN  ANTENNA  HEIGHT: 
STD.  DEV.i 

dev.  uf  mean: 


29.06  METERS 
10.05  METERS 
2 »6G  METERS 


1 : i 


Ti.-ROUT:  NO  VALUES  OUTSIDE  VOUM  2.00  SIGMA  LIMIT 


longitudes  co*«»Ltrfco.  4 Mart  maims  dropped 
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TABLE  2.  KANA  KEOKI  ANTENNA  HEIGHTS  (in  meters)  [Subtract  17  m for  GEOID ] 
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Table  3 

SUVA  ANTENNA  HE I CUTS 


A 


ITE 

: 

DATE 

HT 

12 

1971 

201-202 

46 . 82 

13 

202-207 

59.59 

14 

207-214 

52.01 

28 

1972 

359-360 

58.23 

29 

360-362 

50.29 

All 

Height  from 

above : 

54.33 

3 a ITERATIONS  NO . 


CT 

Zm. 

OBSERV.  HT 

18.58 

8.31 

4 

46 . 82 

16 . 50 

3.  12 

28 

55.66 

17.74 

2.88 

38 

48.87 

20 . 13 

7 . 12 

8 

52 . 32 

16 . 38 

4.54 

13 

50 .29 

17.55 

1.83 

92 

52.  10 

2a  ITERATIONS 

NO. 

a 

Om 

OBSERV 

18.58 

8.31 

4 

14. 11 

2.88 

24 

7.34 

1.64 

20 

12  . 11 

4.58 

7 

16 .38 

4.54 

13 

15.17 

1.63 

87 

B 


12 

45.80 

13 

59.88 

14 

50 . 80 

28 

49.77 

29 

54.59 

All  5 summed: 

54.34 

ITERATION  on 
symmetric  Dopier  cts. 


1 

none 

15 . 88 

3.46 

21 

57 . 04 

17 . 38 

3.41 

26 

51.06 

14.24 

7 . 12 

4 

46.69 

13.86 

4.62 

9 

52.  16 

61 

53.31 

>8  symmetric  Doppler  cts. 


14.99 

3.53 

18 

20.54 

4.96 

17 

15.73 

9.08 

3 

12.60 

4.46 

8 

46 

3a  iterations  by  Satellite 


SAT  -7 

42 

54 

63 

64 

65 


Site 

#13 

52.78 

18.25 

7.45 

6 

72.54 

19.59 

8.76 

5 

66.44 

9.80 

3.70 

7 

56.75 

13 . 70 

5.18 

7 

42.27 

4.94 

2 . 85 

3 

Site  #14 


*55.64 

13.27 

5.02 

7 

63  . 84 

20.43 

8.34 

6 

58.51 

18.21 

6 .07 

9 

53.41 

9.26 

3.50 

7 

33.70 

11.06 

3.70 

9 

* one  value  dropped 


f 


TABLE  4A.  Comparison  Honolulu  Antenna  Heights  (in  meters) 


I 


3a  ITERATIONS 


2a  ITERATIONS 


NO.  OF  NO.  OF 


SITE 

DATE 

HT 

a 

CTm 

OBSERV. 

HT 

a 

am 

OBSERV. 

25 

1972 

306-313 

34.26 

14.54 

3.25 

20 

29.86 

10.05 

2.60 

15 

26 

309-311 

25.62 

13.88 

3.71 

14 

23.91 

12.83 

3.56 

13 

38 

1973 

204-207 

33.16 

25.14 

7.95 

10 

22.66 

10.78 

3.81 

8 

39A 

251-252 

23.04 

26.67 

8.44 

10 

12.27 

11.97 

4.23 

8 

39B 

252-253 

35.08 

27.92 

10.55 

7 

35.08 

27.92 

10.55 

7 

40 

290-297 

18.96 

34.31 

5.11 

45 

19.65 

18.39 

3.20 

33 

41 

1973/74 

353-007 

21.76 

16.81 

1.55 

117 

20.96 

8.94 

1.11 

65 

42 

1974 

9-9 

11.56 

49.11 

24.55 

4 

11.56 

49.11 

24.55 

4 

All  heights  from  above  data  together:  Iterations  on  height  only  (on  HP45)  total  double 

pass  248. 

24.44  18.66  1.23  229  22.58  13.57  0.96  198 


f TABLE  4B.  Honolulu  Antenna  Heights  (in  meters) 

Different  selection  criteria 

For  all  satellites:  elevation  angle  5 to  75°  and  fewer  than  6 

iterations  in  computing  a single  pass  fix 

» 


SELECT 

ION 

SITE 

40 

SITE 

41 

ST.  DEV 

SYM. 

NO.  OF 

NO.  OF 

CTS . • 

HT 

a 

OBSERV. 

HT 

a 

OBSERV 

3a 

all 

18.96 

34.31 

5.11 

45 

21.76 

16.81 

1.55 

117 

2a 

all 

19.65 

18.39 

3.20 

33 

20.96 

8.94 

1.11 

65 

3a 

> 4 

14.31 

17.71 

4.43 

16 

19.94 

15.51 

1.89 

67 

3a 

> 5 

15.49 

17.66 

4.56 

15 

-- 

3a 

> 6 

14.31 

17.75 

4.75 

14 

20.49 

16.02 

2.16 

55 

3a 

> 8 

26.16 

33.51 

9.67 

12 

17.28 

14.42 

2.31 

39 

*(17.82) 

(17.78) 

( 5.36) 

(ID 

3a 

>10 

25.87 

40.27 

14.24 

8 

13.75 

12.25 

2.89 

18 

*(12.71) 

(16.64) 

( 6.29) 

( 7) 

*(15.82) 

( 8.79) 

(2.13) 

( 17) 

*0ne  extreme  value  on  height  dropped  from  line  above. 


TABLE  4C.  Antenna  Heights,  HNL  Site  41  by  Individual  Satellites 


3 ITERATIONS 


ITERATIONS 


NO.  OF 


NO.  OF 


SAT.  NO. 

HT 

a 

CTm 

OBSERV. 

HT 

cr 

CTm 

OBSERV 

42 

23.04 

12.99 

2.71 

23 

21.57 

10.09 

2.45 

17 

34 

23.93 

14.75 

3.08 

23 

19.62 

6.81 

1.60 

18 

63 

-36.11* 

132.55 

50.1 

7 

*same 

64 

29.94 

10.92 

2.33 

22 

25.90 

7.89 

1.91 

17 

65 

13.60 

23.46 

4.6 

26 

10.47 

22.80 

5.53 

17 

99 

17.63 

10.76 

2.35 

21 

13.70 

9.72 

2.51 

15 

FOR  COMPARISON : 

All 

21.76 

16.81 

1.55 

117 

20.96 

8.94 

1.11 

65 

*Two  bad  values  retained,  dropping  the  two: 


39.36 


21.51 


9.62 


t 


TABLE  5.  Antenna  Heights,  Selection  by  Doppler  Counts  and  3a  Iterations 


> 5 Symmetric 
Doppler  Cts. 


> 8 Symmetric 
Doppler  Cts. 


NO.  OF 


NO.  OF 


SITE 

DATE 

HT 

a 

OBSERV. 

HT 

a 

OBSERV. 

CALLAO 

16 

1972 

15-22 

44.6 

33.7 

8.7 

15 

47.4 

37.9 

12.0 

10 

17 

1972 

54-58 

33.2 

31.0 

7.5 

17 

34.1 

25.1 

7.5 

12 

43 

1974 

72-90 

108.7 

53.8 

8.11 

44 

121.1 

75.4 

19.5 

15 

45 

1974 

99-101 

42.6 

54.3 

31.3 

3 

23.3 

-- 

-- 

1 

PALAU 

11 

1971 

167-171 

65.2 

23.0 

6.9 

11 

64.3 

15.0 

6.1 

6 

WELLINGTON 

15 

1971 

241-244 

44.6 

8.7 

2.4 

13 

44.6 

8.7 

2.4 

13 

QUAYAQUIL 

20 

1972 

104-111 

11.1 

20.7 

4.8 

19 

12.0 

23.9 

8.5 

8 

32 

PANAMA 

1973 

50-51 

39.2 

*(15.0) 

55.1 

(12.4) 

24.7 

(6.2) 

5 

(4) 

135.9 

1 

47 

1974 

168-172 

1.3 

*(-4.9) 

28.3 

(16.6) 

7.6 

(4.6) 

14 

(13) 

- 1.4 

31.2 

9.4 

11 

*One  "bad"  double  pass  removed  from  line  above. 


TABLE  6 


SEA  LEVEL- -GEOID  HEIGHT  COMPARISON 


SITE 

GEOID  HEIGHT 

(in  meters) 
SEA 

LEVEL  W. 

RESP.  TO 

SATELLITE  REF. 

SAT  NAV 
MANUAL 

GEM 

6 

SEA  957. 

LEVEL  CONFIDENCE 

<7 

ctM 

No.  of 
OBSERV. 

ANT  HT 
above 
Sea  Level 

GLOMAR  CHALLENGER 

354  DSDP 

-12  to  -15 

-35 

- 42.6 

12.0 

17.9 

5.4 

11 

21.9 

355 

0 to  - 5 

-10 

- 27.4 

30.5 

36.6 

12.9 

8 

21.9 

356 

0 

- 5 

- 4.2 

16.6 

18.1 

6.8 

7 

21.9 

357 

- 2 to  - 5 

-5  to  -6 

- 7.7 

16.4 

13.9 

5.9 

5 

21.9 

358 

- 9 to  -10 

-5  to  -6 

- 10.7 

5.0 

5.9 

2.1 

8 

21.9 

407 

+55 

+55 

50.7 

7.8 

9.4 

3.3 

8 

21.7 

410 

+51 

55  to  60 

53.9 

4.0 

7.3 

1.9 

15 

21.9 

412  (A*) 

42  ± 2 

40  ± 2 

39.80 

16.5 

17.85 

6.75 

7 

22.0 

KANA  KEOKI 

all  Honolulu 

+ 4 to  +5 

0 to  +1 

+ 4.6 

1.9 

13.6 

1.0 

198 

18 

all  Suva 

45  ± 3 

52  ± 2 

34.1 

3.2 

15.2 

1.6 

87 

18 

10  Ponape 

50 

42 

28.0 

24.3 

31.7 

10.6 

9 

18 

11  Palau 

68  to  69 

65  to  66 

51.0 

11.8 

20.4 

5.4 

14 

18 

15  Wellington 

15  ± 2 

12  to  13 

23.7 

5.5 

8.7 

2.5 

12 

18 

20  Quayaquil 

— 

+10 

- 11.9 

6.8 

19.2 

3.3 

33 

18 

34  Antofagasta 

19  to  20 

+30 

39.1 

14.0 

25.2 

6.5 

15 

18 

46  Valparaiso 

20 

19 

36.8 

19.7 

25.6 

8.5 

9 

18 

47  Panama 

0 to  -2 

+5  to  +6 

- 16.1 

7.7 

14.4 

3.6 

16 

18 

*Only  first  part  of  data  used  with  initialized  antenna  height  of  52  meters 
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APPENDIX 


Description 

The  Appendix  contains  all  computer  output  plots  to  support  the  data 
presented  in  Tables 2 through  6.  For  each  site  and  selection  criterion  two 
plots  are  generated.  The  first  plot  shows  all  double  pass  determinations, 
using  selection  criteria  based  on  the  satellite  navigator  output  such  as 
elevation  angle,  number  of  symmetric  Ibppler  count  intervals  and  a particu- 
lar satellite  or  number  of  iterations  required  by  the  satellite  navigator 
to  obtain  the  single  pass  fix. 

The  second  plot  shows  the  remaining  double  pass  determinations  after 
iterations  that  reject  double  pass  results  outside  a specified  multiple  of 

the  standard  deviation,  usually  3 o or  2o.  If  a third  plot  from  a site  is 

I I 

included,  the  no  iteration  limit  is  narrower  (2a)  than  on  the  second  plot 
(3o).  These  iteration  limits  are  indicated  on  the  bottom  line  of  the  plot. 

If  a ()  appears  with  the  iteration  limit,  no  double  passes  have  been 
disregarded  and  therefore  the  plot  is  identical  to  the  first  one.  In  such 
a case  the  first  one  is  not  presented. 

With  each  two  plots  generated  comes  a computer  output  print.  Only  one 
example  of  such  an  output  print  is  given  in  Table  1A  through  1 D for  a 2o 
iteration.  For  more  details  refer  to  the  text. 

The  bottom  line  of  each  plot  contains  the  site  identification,  initialized 
antenna  height  and,  if  applicable,  the  iteration  limits  and  minimum  number 
of  symmetric  Ibppler  counts. 


t 
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Ceneral  order  of  plots  ( identification  by  bottom  line)  : 


Group  I : 


Group  II  : 


Group  III  : 


includes  all  Kana  Keoki  (Pacific  harbor)  sites,  in  numerical 
order  by  increasing  site  numbers,  where  all  double  passes 
with  elevation  angles  between  5°  and  75°  and  up  to  five 
iterations  on  the  satellite  navigator  and  no  restrictions 
on  the  actual  number  or  the  symmetry  of  the  Doppler  counts 
have  been  imposed.  Ifejection  of  double  pass  results  on  a 
3o  and  2o  basis  was  carried  out. 

follows  Group  I and  includes  the  DSDP  sites  (open  Atlantic 
Ocean) . The  output  plots  are  ordered  by  increasing  site 
numbers.  Selection  criteria  are  the  same  as  in  Group  I. 
Eejection  was  carried  out  on  a 2o  basis  only. 

This  group  has  selected  sites  from  Group  I to  show  the 
influence  of  selection  for  increasing  number  of  symmetric 
Ebppler  count  intervals.  The  other  selection  criteria  are 
as  in  Group  I and  Group  II  (refer  to  text  and  Table  3B,  4B, 
and  5).  The  sequence  of  plots  by  site,  selection  and 
rejection  criteria  is  as  follows.  Honolulu  sites  25,  26,  38, 
39A  and  B,  for  six  or  more  symmetric  Doppler  counts  (D.C.>5), 
3a  iterations  : sites  40  and  41,  increasing  number  of  Ibppler 
counts,  3a  iterations. 

Suva  (12,  13,  14,  28,  29) 

Callao  (16,  17,  (A),  (D)) 

Palau  (11) 

Wellington  (15) 

Quayaquil  (20  and  32) 

Panama  (47) 


1 


-A3- 


(all  with  six  or  more  (D.C.  > 5)  and  nine  or  more  ( D.  C.  > 8) 
symmetric  Doppler  counts  and  3o  iterations) 

Group  IV:  This  last  group  shows  the  influence  of  different  satellites  on 

position  and  height  coordinates  (refer  to  text  and  Table  3C  and 
4C).  The  sequence  of  plots  is  : Honolulu  (/Ml),  by  increasing 

satellite  identification  numbers  (which  appear  in  the  plots 
above  the  line  connecting  the  two  x,y  positions  of  a double  pass) 
and  Suva  (#14) , by  increasing  satellite  identification  numbers. 
Selection  criteria  are  as  in  Group  I.  Rejection  by  3o  and  2 a 
iterations  for  Honolulu  and  3a  for  Suva. 

The  results  for  the  computations  are  printed  in  the  three  top  lines  for 
the  mean,  standard  deviation  (a),  and  the  standard  deviation  of  the  mean  (om) 
of  latitude,  longitude,  and  antenna  height. 

Longitude  is  E or  W as  indicated  on  the  right-hand  side  of 
the  X axis. 

Latitude  is  N-positive  or  S-negative. 

Antenna  height  is  calculated  from  double  pass  correction 
added  to  the  initialized  height  (the  latter  is  indicated 
on  the  bottom  line  of  each  plot). 
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IDT..  -1?°  3:7173  •/-  u:0*37  l 81.10  !1  ) OLV.  Ot  flUNi  UTUC03  l 17..'?  fl  J LOT.i  -33°  ?:UI56  •/-  0:0350  l 64. «!  1 i UtV.  Bf  nt'SNi  0:0075  t 13. C3  *1 

L3\S..  77 0 S:37tS  o:ffi*7  l 51.  IB  H ) 0T0P75  t 13.57  M J LONG..  71°  37:fir/0  »/-  Or  05  <7  ( BA.  17  H j 0:0163  ( 75.  Jfc  *1 

nxT.  **lGhT;*5.05  ii  •/-  47. is  n nrv.  a rtiMi  ii.si  n roil.  itlChTi77.su  n «/-  i30.<6  n crv.  oi  ktini  41,75  n 

Trm  °-jaM  tfck) 


■0-Ja’  1 DRY  99-101  1974  KK  CRLLRQ  SITE  »D*  HT:  *38  H # 45  &<r)  DRY  142-146  1974  KK  VALPARAISO  HT:  *39  M # 46  (s^J 


’■3“4  1 DRY  H2-146  1974  KK  VALPARAISO  HT:  *39  M » 46  4a  DRY  .168-172  1974  KK  PRN.  CflNRL  HT: 


•/-  o:o-<3e  ( ei.oe  *i  ) crv.  w meom!  0:0000  ( 12.02  n 3 LnI  , b°  w.is^o  ♦/-  o:o3*e  i n ) dev.  o<  hcrni  0:000?  t n.«i  n 

v.x'j  ♦/-  o:ooou  i ie.?o  n 3 0:0022  ( -t.u7  n 3 long..  ,90  jv.oom  «/-  o:coc3 c n.oo  « 3 0:0016  c 2.eo  n 


DRY  158-172  1974  K K PRN.  CRNRL  HT:  +19  M # 47  3 CT  '°-oc:M  1 qhy  168-172  1974  KK  PRN.  CRNRL  HT:  +19  M #47  ^(f 


l0o511D8fi0  ♦/-  01233S  ( 132.fi.?  .1  3 DEV.  OF  MEBNi  0:0739  ( 135.87  n J 
LOMS..  «"  212521  •/-  01 2251  ( 110.17  M 1 OllOOe  ( 163.13  II  1 

RNT.  HEIGHT*  337.31  n */*  517.11  K OEV.  OF  HERN*  211.81  8 

a-XJ4  T •<  3 


S‘*J<  1 DRY  295-296  DSDP  SITE  353  HT:  -15  M 


•o.xi* 


r 227-23 


-o.ooie 


lst.:  21°  ie:eo20  */-  c E2.R7  n j rev.  re  refiN:  0:0074  c i3.78  h 3 lot.;  21°  ib:796<  •/-  0:0217 1 40.17  n ) rev.  re  reftt  o:oo3E  t 7.10  n 

tree..  157052:0233  ■>/-  0:0215  t 37.19  n 3 oiocsi  i e.77  n j tree..  157052:0239  •/-  o:x77  c 13.3*  n 1 0:0019  c 3.33 1 

***6HT* 2,'£3  11  M-13  11  SEV-  cr  "1E,W‘  ■,-i2  n bnt.  (cigmTi i4.ji  n .7-  i7.7i  n rev.  re'  rermi  «.«3  >i 

f-FCK)  o.x24  TpCK) 


"0'Ja"  1 DRY  290-297  1973  KK  HNL  HT:  *24  M # 40  ^ D.CT^ 


1 CRY  290-297  1973  KK  HNL  HT:  *24  M #40  DOS'  DRY  290-297  1973  KK  HNL  HT:  *24  M #40  ' 1<S- 
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09Y  290-297  3 973  KK  HNL  HT:  *24  H » 40  $ (3ff)  DAY  290-297  1973  KK  HNL  HT:  *24  H # 40  T>.C.>tO  (Jr) 
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3 KK  HM.  HT:  *19  H # 41  D e vi 


DRY  353-372  1973  KK  HNL  HT:  *19  M #41  T).C>U  3cT 
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or  ntw,  0:00*3  ( 9.03  » ) 
0:0033  i 3.v>  n 

*.30  n 


DRY  207-214  1971  KK  SUVR  HT:  *54  M #14  ‘D.C.^  fee) 


vr> 


B l\  5 8 B fi  3 


B 11 

? ?■ 


B 


12"  3:2313  */-  (KI021  l ICO. Si  fl  1 OEV.  OF  ICfW.  0:0107  ( 19.79  M ) LOT..  -12<»  3:J030  ♦/-  0:0922  ( 99.22  tl  ) DEV.  OF  tun,  0:0096  ( >0.31  n 

7°  S' 3733  •/-  KlSSO  t 2fcS.3S  0 > 0:0233  l 42.22  H ) LONG. . /7*>  »:99E7  ♦/-  0:0223  ( 40.33  M > 0:0C3<  t G.Oe  •>  > 


'xsa’  1 DRY  99-101  1974  KK  CRLIRO  SITE  ’D’  HT:  *38  M » 45  '°‘0C*  DRY  99-101  1974  KK  CRLLRO  SITE  *D’  HT: 


LBT.i  /«  19:8144  ♦/-  0:0143  ( 26. *3  *1  ) DEV.  OF  MEAN:  0:0041  ( 7.6S  (1  J 

UST.i  /•  13:8243  -/-  OiOieS  ( 34.47  (1  3 OEV.  OF  MEfNi  0:0040  t 7.35  M 3 LONG..  |34°  27:3967  */-  0:0033  t 17.1!  n 1 <3:0036  t 6.93  B 3 

L»C..  13'"  27:3382  ♦/-  0:0125  t 22.92  M 3 0:0036  t 6.91  M 3 rut.  HElGHTiB4.28  >1  ♦/-  14.37  (1  OEV.  OF  OEflSi  6.11  n 

6NT.  HEIGHT  j 65. 13  fl  ♦/-  22.98  *1  OEV.  01  MEIlNi  6.93  M O-022''  T PC  K ] 

*■**•  tFCK} 


| 
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DRY  167-171  1971  KK  PALAU  HT:  *78  n #11  DOS' 


1973  KK  GUDYHOUIL  HT:  *28  M #32  OC'8  (3»-) 


DRY  168-172  1374  KK  PRN.  CRNRL  HT:  *19  M # 47  DCzZ 
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